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Abstract. The BH mass (and the related Eddington ratio, I = Li,oi / Ledd) in broad line AGN is usually evaluated by 
combining estimates (often indirect) of the BLR radius and of the FWHM of the broad lines, under the assumption 
fT^ ' that the BLR clouds are in Keplerian motion around the BH. Such an evaluation depends on the geometry of the 

BLR. There are two major options for the BLR configuration: spherically symmetric or "flattened" . In the latter 
case the inclination to the line of sight becomes a relevant parameter. This paper is devoted to evaluate the bias 
on the estimate of the Eddington ratio when a spherical geometry is assumed (more generally when inclination 
effects are ignored), while the actual conflguration is "flattened", as some evidence suggests. This is done as 
a function of luminosity and redshift, on the basis of recent results which show the existence of a correlation 
between the fraction of obscured AGN and these two parameters up to at least z=2.5 (date at larger redshifts 
, being insufficient.) The assumed BLR velocity field is akin to the "generalized thick disk" proposed by Collin et 

' al. (2006). Assuming an isotropic orientation in the sky, the mean value of the bias is calculated as a function of 

\^ \ luminosity and redshift. It is demonstrated that, on average, the Eddington ratio obtained assuming a spherical 

. geometry is underestimated for high luminosities, and overestimated for low luminosities. This bias converges for 

all luminosities at z about 2.7, while nothing can be said on this bias at larger redshifts due to the lack of data. 
The effects of the bias, averaged over the luminosity function of broad line AGN, have been calculated. The results 
imply that the bias associated with the a-sphericity of the BLR make even worse the discrepancy between the 
O ' observations and the predictions of evolutionary models. 

^ ' Key words. Galaxies: active - galaxies: nuclei - quasars: general 
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1. Introduction procedure is regarded to be rather safe, but in fact there 
; ^ are still uncertainties on the luminosity and/or redshift 
^ The accretion rate and the black hole mass are dependence of the bolometric correction, 
the two fundamental parameters in our understand- 
ing of the Active Galactic Nuclei (AGN) phenomenon. For the BH mass, two "direct" methods have been fol- 
Measurements of these two quantities are, unfortunately, lowed. The first, applicable only to AGN, is the reverber- 
not devoid of significant uncertainties. ation mapping (RM) method. This method is based on 
The accretion rate m is derived from the bolometric the principle (Blandford & McKec 1982) that the delay 
luminosity, L^ou under assumptions on the efficiency r] in the response of the lines from the Broad Line Region 
for the conversion of gravitational energy, Lboi = -qmc^ . (BLR) to variations of the continuum is a measure of the 
For non-rotating black holes (BH), 77 = 0.057 is generally size of this region, Rblr- Assuming that the Hue widths 
adopted, assuming that effective (for the observer of the are due to motions governed by the BH, the combina- 
electromagnetic radiation) conversion takes place down to tion of Rblr and a velocity derived from the line profiles 
the marginally stable circular orbit at three times the yields a "Keplerian" estimate of the BH mass. The other 
Schwarzschild radius, 37?^; for rotating BH, 77 can reach method is based on the fairly strict correlation between 
the maximum value of 0.42. Lfcoi is generally obtained from the mass of the BH and properties of the stellar bulge 
the luminosity observed in a given band, muhiplied by a of the host galaxy. This method, having proved very reli- 
factor based on the Spectral Energy Distribution (SED) able for a rather large sample of galaxies (Gebhardt et al. 
attributed to the specific class the AGN belongs to. This 2000a, Ferrarese & Merritt 2000, Tremaine et al. 2002), 
represents a benchmark for the previous method when it 
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et al. 2001, Onken et al. 2004). The agreement found, al- 
though far from perfect (e. g. Colhn et al. 2006), has en- 
couraged the extension to many more AGN (especially the 
distant and more luminous ones), for which both the RM 
and the bulge methods can hardly be applied, of a "sec- 
ondary" method. The latter is based on the estimate of 
Rblr through an empirical correlation between this quan- 
tity and the luminosity (see Sect. 2) which has emerged 
from the RM measurements. 



The BH mass can be univocally converted into the 
Eddington luminosity, L^dd- Thus a quantity I can be de- 
fined: I = Liyoi/L^dd, which, altough it tells us nothing pre- 
cise about the accretion rate, is of high interest because 
Lgdd is a very significant physical limit. This quantity is 
often referred to as the Eddington ratio. 

Several papers have been devoted to explore the be- 
haviour of ^, in particular as a function of L and the cos- 
mological epoch, namely the redshift z. Among the un- 
certainties and the selection effects which may plague the 
results, the present paper is devoted to point out and eval- 
uate a particular bias, linked to the possibility that the 
spatial distribution of the BLR clouds is far from spher- 
ical, a situation supported by various lines of evidence. 
The evaluation is based on a recent result on the fraction 
of AGN which are photoelectrically absorbed in the X-rays 
(column density Nh > 10^^ H atoms/cm^, and Compton 
thin), which can be summarized as follows. Calling ^ the 
ratio of the absorbed ones to the total, it turns out that ^ 
is a function of Lx (hence of Lt,oi) (Ueda et al. 2003, La 
Franca et al. 2005) as well as of z (La Franca et al. 2005). 
Qualitatively speaking, in the local Universe this frac- 
tion decreases with increasing luminosity; as the redshift 
grows, the anticorrelation remains but it becomes progres- 
sively shallower. If this behaviour is associated with a lu- 
minosity and redshift dependence of the opening angle of 
the absorbing matter, within which the BLR can be ob- 
served, it should introduce a bias on the estimate of the 
BH mass of broad line AGN (AGN 1 for short), when this 
is performed using the RM method and its "secondary" 
extrapolation. To this effect it is important to stress (see 
Fiore et al. 2003, Perola et al. 2004) that the value Nh = 
10^^ H atoms/cm^ works as a good (the exceptions are a 
minority) discriminant between AGN which are optically 
classified as type 1 and type 2. 

The plan of the paper is as follows. In Sect. 2 the re- 
sults on / from the literature are summarized. In Sect. 
3 two lines of evidence, again from the literature, which 
favour a non-spherical distribution of the BLR clouds are 
briefly described. In Sect. 4 the bias associated with the 
a-sphericity of the BLR, and its dependence on Lboi and z, 
as it can be predicted on the basis of the abovementioncd 
finding, is quantified. A discussion follows in Sect. 6. 



2. Estimates of / = Lboi/Ledd- a summary 

If the BLR clouds are orbiting aroimd the BH, the mass 
of the latter can be estimated as 

RblrYblr. (-^-^ 
Mbh = (ij 

where Rblr is the radius of the region and Vblr. the 
velocity of the clouds. From the application of the RM 
method (Blandford & McKee 1982, Peterson 1993) to a 
sizeable number of objects, a correlation between Rblr 
(as estimated from the "delay" in the response of the 
line chosen for this purpose) and the luminosity has been 
found: 

Rblr « (2) 

where a ~ 0.5, and depends somewhat on the band where 
the luminosity is measured, possibly also on which emis- 
sion line is used (Wandel, Peterson & Malkan 1999, Kaspi 
et al. 2000, 2005, Bentz et al. 2006). Relationship Q is 
applied to AGN samples which include objects without a 
RM estimate. 

The most generally used approach is to estimate Vb lr 
from the FWHM of the line profiles: 

Vblr — k x Vfwhm, (3) 

where k is a geometrical factor which depends on the 
shape of the orbits and on their inclination (Krolik 2001). 
If the orbits are randomly distributed in a spherically sim- 
metric distribution, then (Netzer 1990): 

Vblr = \/3/2 x Vfwhm, (4) 
and eq. Q becomes 

j.,rsphere _ '^RbLrVp^, ^ M fr\ 

~ 4G 

As a general remark, it must be noted that different 
authors adopt a different value of k, which is then assumed 
to be constant. Such an assumption is strictly valid only if 
the width of the lines, as observed, is independent of the 
inclination of the BLR configuration with respect to the 
line of sight (namely the configuration is spherically sim- 
metric); or, with an "on average" meaning, if in a sample 
of objects with an isotropic distribution in their inclina- 
tion angle, there were no limit angle for the BLR to be 
observable, or at least if this angle were independent of 
both redshift and luminosity. 

With this cautionary remark in mind, the results ob- 
tained by various authors (Woo & Urry 2002, McLure & 
Dunlop 2004, Warner et al. 2004, KoUmeier et al. 2005, 
Vestergaard 2002, 2004) on different samples can be sum- 
marized in a simple statement: when selection effects in 
flux limited samples are taken into account, there is no 
evidence that / might depend on Lhoi or/and on z. It is 
hard to evaluate the extent to which this result might 
be influenced by the uncertainties on the extrapolation of 
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^ to large values of Lboi and of z. In the future it may 
become feasible to validate this extrapolation, either di- 
rectly via the RM tecnique, or indirectly by comparing the 
mass estimate with another estimate, obtained indepen- 
dently using the "bulge" method. The uncertainty on the 
bolometric correction seems less relevant. Although all the 
abovementioned authors apply a luminosity independent 
correction, the luminosity dependence found by Marconi 
et al. (2004) is unlikely to change significantly the results. 
Quantitatively, while it is apparent that / is characterized 
by a fairly large scatter, and values larger than unity are 
also found, its mean value settles in the range 0.1 to 0.3, 
depending also on the value of k adopted. This applies to 
Liyoi up to about 10^^ erg/s, and to z values up to 5. 



3. Evidence for a non-spherical BLR 

Two lines of evidence in favour of a non-spherical (but 
axisymmetric) distribution of the BLR orbits are briefly 
recalled here. 

The first case concerns AGN which are also radio- 
galaxies or quasars. Wills & Browne (1986) found a highly 
significant anticorrelation between the FWHM of broad 
up lines and the radio parameter R defined as the ratio 
between the fiux densities of the core and of the radio 
lobes. In the relativistic beaming models of radiogalaxies 
(Blandford & Rees 1978, Orr & Browne 1982, Hough & 
Readhead 1989) the difference between core dominated 
(large R) and lobe dominated (small R) sources is at- 
tributed to a difference in orientation: the first are those 
viewed close to the beam axis, the second those viewed at 
larger angles. Thus R can be regarded as an indicator of 
the orientation angle of the beam, and if the beam coin- 
cides with the axis of simmetry of the BLR, the anticor- 
relation just mentioned can be read as indicating that the 
orbits of the clouds are predominantly confined to a plane 
perpendicular to that axis. This analysis has been re- 
fined and the conclusion confirmed by Wills & Brotherton 
(1995). 

There is no direct evidence that the same picture ap- 
plies to the radio quiet AGN as well, but several arguments 
have been proposed in favour of a flattened, rather than a 
spherical, BLR (Rudge & Raine 1999, McLure & Dunlop 
2002, Jarvis & McLure 2006). 

Collin et al. (2006) made the second case on the basis 
of some significant discrepancies between BH masses es- 
timated with the RM method and those estimated with 
the "bulge" method in the same objects. They show that 
such discrepancies are most marked for the objects with 
the narrowest lines, which are likely the objects with their 
fiattened BLR seen almost pole-on. This conclusion is not 
substantially affected by the use of the second moment of 
the line profile, instead of its FWHM, adopted by Collin 
et al. (2006). 



4. A bias on / due to a flattened BLR 

Following Collin et al. (2006), a simple parameterization 
(which they refer to as a "generalized thick disk" ) will be 
adopted. 



Vblr — 



Vfwhm 



2(a2 



1^)1/2' 



(6) 



where i is the inclination angle of the disk axis relative 
to the line of sight. When ^ is inserted into eq. 1^, an 
estimate of the BH mass is obtained, called here M'^f^: 



Mdisk 
BH 



RblrVfwhm 
4G(a2 sin^ i) ' 



(7) 



Comparing eq. (jSJ with eq. Q, it is apparent that at 
the angle i,=arcsin — a?) the two estimates give 

the same value (i* ~ 35° for a — 0.1 and ~ 30° for a = 
0.3). It is therefore convenient to introduce the ratio 



Ti ^sphere 



3(a' 



sin^ i) 



(8) 



which is illustrated in Fig. ^ for two values of a, namely 

0. 1 and 0.3. The quantity q decreases rather quickly with 

1, when i < i^,, and is fairly sensitive to the value of a. For 
instance, for i = 10°, q is equal to 0.12 (a = 0.1), or to 
0.36 (a = 0.3). When i > i^,, q increases up to a maximum 
of about 3 for both values of a. 



2 - 




Fig. 1. The ratio q = M 



T sphere 



BH 



BH 



function of the 
inclination angle i, assuming a=0.1 (solid line) and a=0.3 
(dotted line). 



As already emphasized by Collin et al. (2006), for a 
single object without an independent observational con- 
straint on its inclination angle, the estimate of the BH 
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mass using eq.© could be in error by up to a factor of 
100. However, the interest here is on systematic errors on 
large samples. 

If the orientation in the sky of the disk axis is assumed 
isotropic, namely 



diV 



= sm I , 



(9) 



it immediately follows that, if there were no limits to the 
angle at which the BLR can be observed, the mean value 
of q would be 2 (a=0.1) or 2.3 (a=0.3). If a spherical dis- 
tribution is adopted (similarly for a k value of 1), Mbh is 
significantly overestimated, and hence I underestimated. 
This is relevant in itself, but before discussing in this re- 
spect the results reported in Sect. 2, one must deal with 
the fact that in the AGN unification model the type 1 ob- 
jects can be found only within a limiting angle, say zq, as 
pointed out by Collin et al. (2006). In addition this angle, 
according to the findings by Ueda et al. (2003) and La 
Franca et al. (2005), is a function of both luminosity and 
redshift. 

With reference to the results of La Franca et al. (2005) 
on the fraction ^ as a function of Lx and z, obtained 
taking into account the selection effects on the samples 
used (illustrated in their Figure 11), one can immediately 
derive an opening angle zq: 



cosio = ^.{Lx, z) 



(10) 



which discriminates type 1 AGN (where the BLR is visi- 
ble) from type 2 AGN (where it is not). The function ^, in 
La Franca et al. (2005), holds in the ranges 0.25< z <2.75 
and 42. 5< \ogLx <45.5. At redshifts and luminosities out- 
side this ranges, ^ is kept costant and equal to the val- 
ues obtained at the limits of the intervals. The outcome 
is illustrated in Fig. [3 where Lx has been converted to 
Lf,oi using the (luminosity dependent) bolometric correc- 
tion given by Marconi et al. (2004). Each line refers to a 
luminosity interval of ±0.5 dex. 

It must be stressed that La Franca et al. (2005) used 
a simple functional form to describe the quantity ^, which 
does probably represent fairly well the real situation up to 
z around 2. However, it cannot be extrapolated to larger 
values without a sizeable change in slope, which however 
could not be properly estimated with the data available. 
The authors preferred to introduce a sort of "saturation" 
value, which explains why the curves for the four different 
values of L^oi intersects at z = 2.7, with — 41°, which 
corresponds to the "saturation" value of ^ = 75%. 

At this point, assuming that the disk axis and the ob- 
scuring matter axis are coincident, the mean value of 
as a function of luminosity and redshift, can be derived 
using eq. Q: 



ia(Lx,z) 




q sin i di 



Jq 



io{Lx,z) 



sin i di 



(11) 



where io{Lx,z) is given by eq. (|10|l . The results are illus- 
trated in Fig.^lfor the values a = 0.1 and 0.3. 
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Fig. 2. The opening angle jq, inferred from the fraction 
in La Franca et al. (2005), as a function of z for different 
values of the AGN bolometric luminosity. 



5. Discussion 

The curves in Fig.^show that the mean value of the bias, 
<q>, is typically greater than 1 and can be as large as 
about 2. Furthermore, this bias increases with the lumi- 
nosity, converging for all luminosities towards a value of 
order unity at z larger than 2. As noted previously, the es- 
timate of the opening angle iq, hence of <q>, ceases to be 
reliable in this regime of the cosmic time, because the lack 
of data prevented La Franca et al. (2005) from evaluating 
the actual evolution of ^ further back in time. The "satu- 
ration" effect adopted in ^ corresponds to an angle ~ 40°, 
such that, according to eq. (|lll) . <q> is ~ 1 for a=0.3, or 
~ 0.7 for a=0.1. Since the weight of the solid angle within 
this value of ig is relatively modest, the mean value of q 
will result significantly lower than unity only for iq much 
lower than i*. If, for instance, ^ were to reach 90% (that 
is io = 26°) at values of z much greater than 2, than <q> 
would be 0.32 (a=0.1) or 0.56 (a=0.3). 

We also calculated the rms of the q distribution and 
found values of 0.9, 0.6 and 0.2 for a=0.1 and io=90°, 
60°, 30°, respectively (the corresponding values of <q> 
are 2.0, 1.3 and 0.4). Similar values are found for a=0.3. 
The expected scatter in masses is thus reassuringly smaller 
than that found in the Black Hole mass (estimated from 
the reverberation mapping) vs. bulge dispersion velocity 
relationship (Onken et al. 2004). 

Turning to the parameter this quantity is underesti- 
mated when <q> > 1, if Mbh is "measured" assuming a 
spherical distribution of the BLR, the other way round if 
<q> were less than unity. 

A possible way to illustrate the effects of this bias 
consists in calculating its value averaged over the entire 
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Fig. 3. Mean value of < g >, as a function of rcdsliift for different values of the bolometric luminosity; a=0.1 (left) 
and 0=0.3 (right). 



Luminosity Function. La Franca et al. (2005) give both 
the Luminosity Function (LF) and the ^ which best fit the 
data (fit #4 in their Table 2), after taking into account 
the selection effects (in the X-ray and optical bands) for 
the samples used. The LF behaviour with cosmic time fol- 
lows a Luminosity Dependent Density Evolution. One can 
therefore combine the LF and ^ to obtain, for a given z, 
the average bias, <q>L, over the entire luminosity range: 



< q >L-- 



JiogLxi diogLx '^^°SLx Jo q sm i dz 



r 



logLx2 d<I>i(L^,z) 
igLxi dlogLx 



dlogLxC^^'^'^sinzdz 



-,(12) 



where $i is the luminosity function of unabsorbed AGN 
(which depends on ^) and the parameter a is included in 
the function q (eq. 0). This result should, by construc- 
tion, be free from selection effects. For a choice of logLxi 
= 42 and logLx2 = 48, <q> l is illustrated in the left panel 
of Fig. 0] (the corresponding values of logLboi, obtained ap- 
plying the luminosity dependent bolometric correction of 
Marconi et al. (2004), are: logiboa — 43.030 and logLf,oi2 
= 50.937). 

To illustrate the effect arising from the inclusion of the 
numerous AGN in the decade around logLx = 41 (logLfcoj 
= 41.893), in the right panel of Fig. 0] (and only for the 
case a=0.1) <q>L is compared with the result shown in 
the left panel of the same figure. The effect is that, the 
higher is the luminosity, the steeper is the dependence of 
the mean bias <q>L on z. 

A direct application of our results to the samples used 
by the various authors, quoted in Sect. 2, is not straight- 
forward, because this would require proper evaluation of 
their specific observational selection effects, a hard task 
which goes beyond the aims of this paper. However, some 
general remarks can be made. The results shown in Fig. 
pimply that, if I were intrinsically constant with z, and 
if its "mean" value were calculated with the spherical ap- 
proximation (or, more generally, with a constant value of 



K in eq. Q) then one should observe an increase of I with 
z. Since the observational results indicate instead a con- 
stant value, the bias discussed in this paper implies that 
the actual value of I decreases with z. 

In semianalytical models which link the evolution of 
the galaxies in the hierarchical clustering scenario with 
the quasar evolution (e.g. Menci et al. 2003, 2004), the 
black hole accretion is triggered by galaxies encounters. 
In this scenario, at high z the protogalaxies grow rapidly 
by hierarchical merging, meanwhile much cold gas is im- 
ported and also destabilized, so that the black holes are 
fueled at their Eddington rates. At lower z the accretion 
rate of cold gas onto the central black hole diminish due to 
the combined effects of the decrease of the galaxies merg- 
ing and encounter rates and the decrease of the amount of 
galactic cold gas, which was already converted into stars 
or accreted onto the black hole. This model predicts an av- 
erage Eddington ratio dropping from Z ~ 1 at z ~ 2.5 to I 
~0.01 at z ~0 (Menci et al. 2003). Our results imply that 
the bias associated with the a-sphericity of the BLR make 
even worse the discrepancy between the observations and 
the predictions of the models. 
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